APS/123-QED 



Quasiparticle Dynamics and Phonon Softening in FeSe Superconductors 
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Quasiparticle dynamics of FeSe single crystals revealed by dual-color transient reflectivity mea- 
surements (AR/R) provides unprecedented information on Fe-based superconductors. The ampli- 
tude of the fast component in AR/R clearly gives a competing scenario between spin fluctuations 
and superconductivity. Together with the transport measurements, the relaxation time analysis fur- 
ther exhibits anomalous changes at 90 and 230 K. The former manifests a structure phase transition 
as well as the associated phonon softening. The latter suggests a previously overlooked phase tran- 
sition or crossover in FeSe. The electron-phonon coupling constant A is found to be 0.16, identical 
to the value of theoretical calculations. Such a small A demonstrates an unconventional origin of 
superconductivity in FeSe. 

PACS numbers: 74.70.Xa, 78.47.D-, 78.47.J- 



Since the discovery of Fe-based superconductors 
(FeSCs) in 2008 Q, tremendous experimental and theo- 
retical effort has been devoted to exploring their charac- 
teristics. These Fe-based pnictide compounds exhibit a 
very interesting phase diagram, with antiferromagnetism 
(or spin-density wave) at low doping and superconduc- 
tivity at intermediate doping Q . The simultaneous pres- 
ence of magnetism and superconductivity in the phase 
diagram implies that magnetism plays an important role 
in the superconductivity mechanism. The existence of 
precursor superconductivity above T c which competes 
with the spin-density wave order and a pseudogap- 
like feature with onset around 200 K [4| were observed on 
underdoped (Ba, K)Fe2As2 and nearly optimally doped 
SmFeAsOo.8Fo.2, respectively. Additionally, a coherent 
lattice oscillation was also found in Co-doped BaFe2As2 
using time-resolved pump-probe reflectivity with 40 fs 
time resolution |5|. Among various FeSCs, the iron 
chalcogenide FeSe [6j stands out due to its structure sim- 
plicity, which consists of iron-chalcogenide layers stack- 
ing one by another with the same Fe +2 charge state as 
the iron pnictidcs. This so-called "11" system is so sim- 
ple that it could be the key structure to understand- 
ing the origin of high- T c superconductivity Q ■ There 
has been considerable concern over the interplay between 
electronic structure, phonons, magnetism, and supercon- 
ductivity in 11-type FeSe. Therefore, further studies of 
their quasiparticle dynamics are indispensable to under- 
standing the high- T c mechanism in FeSCs. Here we 
report the first time-resolved femtosecond spectroscopy 
study of FeSe single crystals to elucidate the electronic 
structure and the quasiparticle (QP) dynamics. 

In this study, FeSe single crystals were grown in evac- 



uated quartz ampoules using a KCI/AICI3 flux [8j. The 
crystalline structure of the samples was examined by x- 
ray diffraction. The superconducting transition temper- 
ature T c of the FeSe single crystal was determined to be 
8.8 K by the middle point of the resistive transition. The 
femtosecond spectroscopy measurement was performed 
using a dual-color pump-probe system (for light source, 
the repetition rate: 5.2 MHz, the wavelength: 800 nm, 
and the pulse duration: 100 fs) and an avalanche pho- 
todetector with the standard lock-in technique. The po- 
larizations of the pump and probe beams were perpen- 
dicular to each other and parallel to the ah plane of FeSe 
single crystals with no distinction between a and b ori- 
entation. The fluences of the pump beam and the probe 
beam are 9.92 and 1.40 /xJ/cm 2 , respectively. The pump 
pulses have corresponding photon energy (3.1 eV) where 
the higher absorption occurred in the absorption spec- 
trum of FeSe Q and hence can generate electronic exci- 
tations. The QP dynamics is studied by measuring the 
photoinduced transient reflectivity changes (AR/R) of 
the probe beam with photon energy of 1.55 eV. 

Electronic excitations generated by the pump pulses 
result in a swift rise of AR/R at zero time delay. The 
observed excitation is triggered by transferring the elec- 
trons from d valence band of Fc to d conduction band 
of Fe [l0l | . At zero time delay, the number of the ex- 
cited electrons generated by this nonthermal process is 
related to the amplitude of AR/R. These high-energy 
electrons accumulated in the d conduction band of Fe re- 
lease their energy through the emission of longitudinal- 
optical (LO) phonons within several picoseconds (Tll[l2T|. 
The LO phonons further decay into longitudinal-acoustic 
(LA) phonons via anharmonic interactions, i.e., transfer- 
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FIG. 1. (color online) (a) Temperature and delay time depen- 
dence of AR/R in an FeSe single crystal. The inset shows a 
part of the AR/ R on an enlarged scale. Solid lines indicate 
AR/R= 0. (b) Temperature dependence of the resistivity p 
indicates the high quality of an FeSe single crystal. The kink 
of p(T) at 90 K manifests the structure phase transition. 




40 80 120 160 200 240 

Temperature (K) 



' (b) 1 

m 2 ' 
3 '■ 






5 6 7 8 9 10 

Temperature (K) 













8 




A 


S B 

5 






3 







Temperature (K) 



40 80 120 160 200 240 280 

Temperature (K) 



FIG. 2. (color online) Temperature dependence of the ampli- 
tude (a) A c , (c) Alo and the relaxation time (b)r e , (d) tlo 
by fitting Eq. JTJ). Inset of (a): Density of thermally excited 
QPs, m{T). A solid line fits to the Rothwarf- Taylor model. 
Inset of (b): The temperature-dependent r c is shown between 
4 and 10 K. Dashed lines are guides to the eyes. The arrow 
indicates the divergence of r c around T c . 



ring energy to the lattice. This relaxation process can 
be detected using a probe beam as shown in Fig. 1(a), 
which shows the 2D AR/ R taken on an FeSe single crys- 
tal. Four temperature regions appear. Above 230 K, T* 
(region I) , there is a fast negative response with a relax- 
ation time of about 1.5 ps together with a long period 
of oscillation. When the temperature decreases into re- 
gion II, a positive and slow response appears and AR/ R 
gradually becomes smaller until T = 90 K (T s ) where 
the resistivity p(T) curve shows the kink [Fig. 1(b)]. 
Below 90 K (region III), the slow positive response dis- 
appears and is replaced by a complicated mixture of the 
positive and negative components as discussed later. For 
T < T c (region IV), a long-lived negative response ap- 
pears similar to the one in region I. These four regions in 
2D AR/ R are apparently corresponding to the feature of 
p(T) in FeSe. The relaxation processes (t > 0) of AR/R 
in FeSe can be phcnomcnologically described by 

^ = A e e-*' T ' + 4oe- (/ -° + A Q 
K 

+A LA e- t/T ^ A sm[2irt/T(t) + 0] (1) 

The first term in the right-hand side of Eq. (JTJ) is the 
decay of the excited electrons with an initial population 
number, A c , and a relaxation time, r e . The second term 
is the phonon number, A^o, and a corresponding decay 
time, tlo- The third term describes energy loss from 
the hot spot to the ambient environment within the time 



scale of microsecond, which is far longer than the period 
of the measurement (~ 150 ps) and hence is taken as a 
constant. The last term is the chirped oscillation compo- 
nent associated with strain pulse propagation [l3l ]: A la 
is the amplitude of the oscillation; tla is the damping 
time; T(t) is the time-dependent period; (f> is the initial 
phase of the oscillation. 

We first discuss the temperature-dependent AR/R in 
FeSe. Fig. QTa) shows AR/R from 290 down to 4.4 K, 
from which each component described above can be ex- 
tracted using Eq. (JlJ. Results of the extraction are 
shown in Fig. (5Ja)-2(d). For the fast component in 
AR/R, a remarkable increase in the amplitude (^4 ) is 
shown in Fig. Ufa) when T decreases from 90 K. The 
temperature dependence of spin fluctuations below T s 
in FeSe which were revealed by 77 Se NMR measure- 
ments [l5| is totally different from that in cuprates [r|. 
Therefore, this temperature-dependent A c below T s is 
attributed to the strong antifcrromagnctic spin fluctu- 
ations which may associate with the structural trans- 
formation at 90 K rather than with the pseudogap in 
cuprates. Intriguingly as T < T c (region IV in Fig. [T]), 
A e decreases due to the onset of superconductivity. This 
suggests that the spin fluctuations and superconductivity 
are competing factors in the FeSe system. The above re- 
sults provide the first compelling experimental evidence 
of competing orders in FeSe, which are consistent with 
the theoretical calculations [T7[. It is noted that the ex- 
perimental evidence on the competing orders was also 
reported in underdoped (Ba, K)Fe2As2 0- According to 
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the Rothwarf- Taylor model (l8| , the density of thermally 
excited QPs n T cx [A C {T) / A C {T - 1 as shown in the 
inset of Fig. . The temperature-dependent behavior 
of n T can be further fitted by n T oc [A( T) T] 1 / 2 exp[- 
A(T)/T], where A(T) is the superconducting energy 
gap. Assuming a BCS temperature dependent A(T) = 
/\{<S)[l-{T/T c )] 1 / 2 , the fits lead to A(0) = 3.24/c B T c in 
accord with the value obtained from specific heat mea- 
surements [1]. [This isotropic A(T) scenario is certainly 
too simplified, but it suffices to estimate the scale of 
A(T).] 

The presence of a gap in the QP density of states 
gives rise to a bottleneck for carrier relaxation, which 
is clearly observed in the relaxation time r c close to T c . 
The mechanism of the bottleneck can be described by 
the Rothwarf- Taylor model [HI]. When two QPs with 
energies higher than A, a high energy boson (HEB) with 
energy u> > 2A is created. The HEBs that remain in the 
excitation volume subsequently break additional Cooper 
pairs, effectively preventing QP from recombination. Till 
lj < 2A and the Cooper pairs could not be broken further 
by HEBs, the number of QPs finally decreases in several 
picoseconds. In the case of a mean-ficld-likc gap, i.e., 
the gap gradually shrinks with T approaching T c from 
the low-temperature side, more HEBs are available to re- 
generate QPs. The recombination processes become less 
and less efficient. Hence, r below T c through Cooper 
pair recombination is longer than 1.5-2 ps and diverges 
around T c , which is dominated by the superconducting 
gap opening. To compare, r is ~ 1.5-2 ps between T c 
and 90 K and is dominated by the appearance of spin 
fluctuations in this region. 

For the slower component, the amplitude ^4lo 
monochromatically decreases as T does, and then it com- 
pletely disappears in the superconducting state. On the 
contrary, the relaxation time tlo exhibits two marked 
anomalies both at T*~ 230 K and T s ~ 90 K, corre- 
sponding to the boundary between region I and II and 
between region II and III, respectively. During the struc- 
tural transition from a tetragonal phase to an orthorhom- 
bic phase [3] , the energy of LO phonons more efficiently 
releases to lattice and results in a significantly shorter 
relaxation time tlo at TV On the other hand, the ab- 
normally long relaxation time tlo around 230 K sug- 
gests an elusive higher-temperature phase transition or 
crossover to undermine the energy release efficiency. The 
sign change of Seebeck coefficient of FeSe was found to 
be also at T* (20|. 

By fitting the AR/R curves with Eq. (TTJ), dynamic 
information on QPs and phonons is available, which in- 
cludes the number of QPs, the relaxation time of QPs, 
and the energy of phonons. In a metal, the photoinduced 
QP relaxation time is governed by transferring energy 
from electron to phonon with electron-phonon coupling 
strength A 21 1. 
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FIG. 3. (color online) Temperature-dependent peaks and 
depths of the oscillation component in AR/R. Inset: Temper- 
ature and delay time dependence of the 3D oscillation compo- 
nent was obtained by subtracted the decay background [the 
1 st , 2 nd , and 3 rd terms in Eq. (fTJ] from AR/R of Fig. [TJa). 
Solid lines are guides to the eyes. 
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where X(u> 2 ) is the second moment of the Eliashberg func- 
tion and T e can be further described by 22 1 



T P = 




(3) 



where R is the unperturbed reflectivity at 400 nm, F is 
the pumping fluences and 7 is the linear coefficient of 
heat capacity due to the electronic subsystem. The mean 
value is taken for the depth z going from the crystal 
surface down to the skin depth l s ~ 24 nm (estimated 
from the skin depth of electromagnetic wave in metal, 
X/Ank). All of the parameters for the calculations of 
electron-phonon coupling strength were listed in Table 
I. For the estimate of (w 2 ), some vibrational modes 
arc more efficiently coupled to QPs than others are. 
In the case of Co-doped BaFc2As2, the symmetric Ai g 
mode is coherently excited by photoexcitation and 
efficiently coupled [5[. Consequently, we take the Ai g 
mode into account in the present case of FeSe, which is 
the strongest phonon mode in electron-phonon spectral 
function, a 2 F(u) [Io| . By Eq. ©, the consequent 
electron-phonon (Ai g mode) coupling constant is A = 
0.16 in FeSe. This value is consistent with the theoret- 
ical results of A = 0.17 [l(| obtained by using linear 
response within the generalized gradient approximation 
(GGA). Furthermore, we can use the McMillan formula 
T c =((w)/I.2)cxp{-[1.40(1 + A)]/[A - (x*(l + 0.62A)]} 
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TABLE I. The parameters for A estimated at T= 


20 K of an 


FeSe single crystal. 




Tc (K) 


R (400 nm) F (fij/cm 2 ) 7 a (mJmor'K" 2 ) 


T e (ps) 


A b lg (meV) X(co 2 ) c (meV 2 ) 


A 


FeSe 8.8 


0.25 9.92 5.73 


1.75 


19.9 61.3 


0.16 


a From Ref. @ 


. b From Ref. [H]. c Obtained from Eq. ©. 









to evaluate the critical temperature T c [23[. Taking 
(hkj)= 19.9 meV and 0, we obtain T c ~ 0.08 K, 
which is one order of magnitude lower than the actual 
T c of about 8.8 K. Actually, the total coupling strength 
between quasiparticles and bosons is estimated to be 
1.55 from the specific heat measurements [1]. Therefore, 
the present study experimentally verifies that the 
superconductivity in FeSe should be attributed to the 
unconventional mechanism other than electron-phonon 
coupling. 

Further insight into the phase transition in FeSe is pro- 
vided by the study of the oscillation component of AR/ R. 
Temperature dependence of a strain pulse (LA phonons) 
propagation was clearly observed in the oscillation fea- 
ture of AR/R after subtracting the decay background 
[i.e., the 1 st , 2 nd , and 3 rd terms in Eq. Q], as shown 
in the inset of Fig. [3l In order to quantify the proper- 
ties of LA phonons in FeSe, the locations of the peaks 
and depths of the oscillation in the inset of Fig. [3] are 
shown as the function of temperatures in Fig. [3] At 
high temperatures, the damping time is very short and 
the oscillation sustains only for one period. However, the 
number of oscillation periods significantly increases be- 
low 100 K; hence the damping time becomes much longer. 
This means that the LA phonons can propagate further 
into the interior of FeSe crystals with the orthorhom- 
bic structure. According to the difference between 1 st 
depth (at 23.52 ps) and 2 nd depth (at 72.78 ps) at T= 
110 K, the phonon frequency is found to be 20.3 GHz. 
The phonon energy is estimated to be ~ 0.087 meV. The 
coherent acoustic phonon detected by a pump-probe re- 
flectivity measurement can be described as a Brillouin 
scattering (24T | phenomenon occurring in the materials 
after the excitation of pump pulses. The scattering con- 
dition is Qphonon = 2nfc prob o cos(6>i) , where g p honon is the 
phonon wave vector, n is the real part of the refractive 
index, and the probe photon has a wave vector fc pro be 
arriving at an incident angle 9{ (inside crystals) with 
respect to the surface normal. Following this scatter- 
ing condition, the probe beam acts as a filter to select 
the acoustic wave propagating along the scattering plane 
symmetry axis, i.e., the normal to the crystal surface and 
traveling with the wave vector g p honon- The energy of 

the aCOUStic wave is iJphonon = ^phonon = ft<7phononW S = 

h2nv s k pro he cos((9i), where v s is the sound velocity along 
normal direction of crystal surface. Using A pro bc= 800 
nm, n pro b c = 2 Q, 0\— 2.5 ° [estimated from the incident 
angle (5°) of the probe beam by Snell's law] and v s = 3.58 
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FIG. 4. (color online) Temperature dependence of the phonon 
energy derived from the oscillation component in Fig. [3] The 
inset shows a part of the temperature-dependent phonon en- 
ergy on an enlarged scale. Dashed lines are guides to the 
eyes. 



km/s [25|, the phonon energy, Sphonom is calculated to be 
0.077 meV, which is very close to the result, 0.087 meV, 
directly obtained from the above AR/R measurements. 
We further investigate the temperature dependence 

of the LA phonon energy as shown in Fig. HJ The 

phonon energy dramatically drops by 60% |26| around 
90 K where a structural phase transition occurs, and 
then keeps constant at low temperatures. Additionally, 
we found that the LA phonons also soften by 6% in the 
superconducting state as shown in the inset of Fig. SJ 
which is consistent with the larger distance between l Bt 
depth and 1 st peak in Fig. [3] Very recently, this soft- 
ening of lattice has been found in Co-doped BaFe2As 2 
by resonant ultrasound spectroscopy [27], HI]. Fernan- 
des et al. [27]| found the 16% softening of shear modulus 
in BaFei.84Coo.i6As2 at T c = 22 K. For the nonsupcr- 
conducting case of BaFc2As2, however, the rather large 
softening of 90% was observed around 130 K, which is 
the structural and antiferromagnetism phase transition 
temperature. Similarly, a large phonon softening due to 
structural phase transition and a rather small phonon 
softening due to the superconductive phase transition 
was also observed in 11-type FeSe in this Letter. These 
results suggest that the reduction of phonon energy at 
both the structural and the superconducting phase tran- 
sition is a general feature in FcSCs. The above phonon 
softening may participate in the superconducting pair- 
ing, albeit not the mechanism responsible for high T c in 
FcSCs. 
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In summary, we have studied the ultrafast quasiparti- 
cle dynamics and phonon softening in FeSe single crystals 
by dual-color femtosecond spectroscopy. The tempera- 
ture dependence of the amplitude A suggests a com- 
peting scenario between spin fluctuations and supercon- 
ductivity. The relaxation time r c of AR/R reveals an 
electron-phonon coupling constant A= 0.16. The struc- 
ture phase transition at 90 K is elucidated through both 
the electrical transport properties and the anomalous 
changes of relaxation time tlo of AR/R. A previously 
unknown feature at 230 K is further identified. Moreover, 
the energy of LA phonons at 110 K was estimated to 
be 0.087 meV from the oscillation component of AR/R, 
which markedly softens around both the structural phase 
transition and superconducting transition. Our results 
provide vital understanding of the competing picture be- 
tween the spin fluctuations and superconductivity and of 
the role of phonons in Fe-based superconductors. 
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